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PACS 7 5 . 6 . E j ~ Magnetization curves, hysteresis, Barkhausen and related effects 

Abstract. - We obtain a general formula for the distribution of sizes of "static avalanches", or shocks, in 
generic mean-field glasses with replica-symmetry-breaking saddle points. For the Sherrington-Kirkpatrick 
(SK) spin-glass it yields the density p{AAI) of the sizes of magnetization jumps AAI along the equilibrium 
magnetization curve at zero temperature. Continuous replica-symmetry breaking allows for a power-law 
behavior p{AAI) ~ l/(AAf )^ with exponent t — 1 for SK, related to the criticality (marginal stability) of 
the spin-glass phase. All scales of the ultrametric phase space are implicated in jump events. Similar results 
are obtained for the sizes S of static jumps of pinned elastic systems, or of shocks in Burgers turbulence in 



large dimension. In all cases with a one-step solution, p{S) ~ ^e 



A simple interpretation relating 



droplets to shocks, and a scaling theory for the equilibrium analog of Barkhausen noise in finite-dimensional 
spin glasses are discussed. 
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Many disordered systems crackle when driven slowly, reacting 
with abrupt responses over a broad range of scales |T|. These 
avalanche phenomena occur in granular materials I^J, earth- 
quakes (3), fracture Q, liquid fronts ||5|, vortex lattices |J6|, and 
other pinned elastic objects such as domain walls in disordered 
ferromagnets [71, where jumps in magnetization are known as 
Barkhausen noise ||8l. In electronic glasses, the analogue of a 
magnetic hysteresis experiment is gating, which exhibits strik- 
ing memory effects L9J. Also in that context one expects crack- 
ling upon increasing the carrier density. 

The size S of these events is power-law distributed, i.e. scale- 
free, p{S) ~ S"^^. This property, often termed self-organized 
criticality, emerges naturally in sandpile models, where analyt- 
ical results were obtained [10|. But even there, p{S) is dif- 
ficult to compute. Scale-free response also occurs in pinned 
elastic systems, where quenched disorder leads to glassiness 
and metastability at all scales. The distribution of avalanche 
sizes for a single elastic interface was obtained from the func- 
tional renormalization group (FRG) |il|, and compared with 
numerics p2[[T3] and with wetting experiments at the depin- 
ning transitionpj . The random-field Ising model with short- 
range interactions, much studied in this context, exhibits a tran- 
sition between non-critical and infinite avalanches as disorder 
is varied 1 14 15), with scale-free avalanche distributions only 
at a special point in the phase diagram. While domain-wall 
motion plays an important role in soft magnets, a description 
without nucleation, long-range dipolar interactions and the en- 



suing frustration between domains would be incomplete 1 8 1 . 
The situation is less explored in strongly frustrated spin- 
glasses, whose complex energy landscape shares many features 
with that of pinned elastic systems. In particular, the spin-glass 
phase exhibits criticality with power-law spin correlations, as 
predicted in mean-field theory ||16J and in the droplet picture 
1 17 1. This property is difficult to access by standard experimen- 
tal protocols. However, the statistics of magnetization bursts in 
a hysteresis experiment (the Barkhausen noise) should be sen- 
sitive to the criticality of the glass state, and thus serve as a 
probe of spin-glasses, both experimental and numerical. 
The aim of this Letter is to compute the statistics of equi- 
librium (i.e. static) magnetization jumps in the Sherrington- 
Kirkpatrick (SK) mean-field spin-glass. We obtain a formula, 
(|7| below, which applies more generally to any mean-field 
model described by a replica-symmetry breaking (RSB) sad- 
dle point. The strategy is similar to Ref. 1 1 1 1 for elastic inter- 
faces: A static avalanche, or shock, occurs when the system 
jumps discontinuously between two degenerate global minima 
as the energy landscape is tilted with an external force. This 
phenomenon shows up as non-analytic cusps in all moments 
of the effective (i.e. renormalized) force. Calculation of these 
cusps allowed 1 1 1 1 to obtain the jump-size distribution in an e 
expansion around internal dimension rf = 4. To lowest order, 
it was found to be identical for dynamical avalanches (i.e. un- 
der a driving force) and for static avalanches, or shocks, with 
T — 2 — j^ + O(e^), C being the roughness exponent. Here, 
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we extend previous mean-field studies of the second moment 
of the equilibrium effective force 1 18-20) to all moments. The 
resulting formula (JTl) allows for a simple interpretation. We 
suggest an extension as a scaling theory for finite-dimensional 
spin glasses. It is based on a relation between shock and droplet 
distributions, which extends an identity obtained for a particle 
pTj . Similar shocks occur in many systems, and we consider 
elastic manifolds [22[ as well as decaying Burgers turbulence 
[23) , embedded in large dimension N, the mean-field limit. We 
find that one-step RSB always results in p{S) ^ Se~^^ , sim- 
ilar to the shock-size distribution obtained by Kida in Burgers 
turbulence |J24][25i|. In contrast, we show that continuous RSB 
results in novel and interesting scale-free avalanche distribu- 
tions. 

Mean-field spin glasses and the SK model - Out of equilibrium 
avalanches in the SK model at T = were studied numeri- 
cally p6) , and found to exhibit criticality, i.e. power-law size 
distributions. The system self-organizes to remain always at 
the brink of stability. Upon increasing the external field by 
typically ^ l/\fN, a first spin flips, and with finite probabil- 
ity entails 0{N) spin flips, with a change in magnetization of 
0{y/N). The thermodynamic criticality of the spin-glass phase 
suggests similar avalanche phenomena at equilibrium. Indeed, 
it has long been known 1 27 1 that the total equilibrium magne- 
tization M{h) of the SK model undergoes a sequence of small 
jumps, akin to mesoscopic first-order transitions, as the field h 
is increased. These jumps are of size N-^^^ and lead to non- 
selfaveraging spikes in the susceptibility. However, the ana- 
lytical understanding of avalanches in spin-glasses, and their 
relation to thermodynamics, has remained scarce. 
As a first step |28|, the equilibrium magnetization of mean- 
field systems with p-spin interactions (p > 2) was analyzed 
and compared to a toy model of a large set of states with ran- 
dom energies E and magnetizations M |J29]|. When the free 
energies Ei,2 — hMi,2 of the two lowest states cross as h is in- 
creased, ajump |Mi — A/2I in magnetization occurs. This basic 
picture, with metastable states replacing the random states, re- 
mains a good qualitative guide, both in the p-spin model and in 
the SK model of spin glasses, as we see below. In Ref. flE], 
the presence of large avalanches in the equilibrium magnetiza- 
tion was proved by exhibiting a non-analyticity in the second 
moment [Af (/12) - M(/ii)]2 for /12 - /ii ~ N''^/'^, in close 
analogy to the force correlator of elastic systems pO) . Such 
shocks are sharply defined only at T = 0. Thermal smearing of 
the magnetization curve buries the presence of rounded jumps 
as soon as T > l/\/]V, as the calculation confirms. Thus, in 
mean-field spin models, including SK, very low temperatures 
must be considered. 
Consider now specificaUy the SK model, described by (p = 2): 
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where the J^j are i.i.d. centered Gaussian random variables of 
variance J"^ /N, that couple all TV Ising spins, and h is the ex- 
ternal field which will be varied adiabatically. This problem is 
more involved than the p-spin model with p > 2, since its glass 



phase involves infinite-step RSB with marginally stable states, 
unlike p > 2, which has a 1-step solution. It is reflected in 
crucial differences in the avalanche statistics. The equilibrium 
solution of ([T]l at A^ — > 00 is given by Parisi's full replica- 
symmetry breaking ansatz for the saddle point of the overlap 
matrix Qab ~ {<7a<^b)- The order parameter is a monotonous 
function q(x) on the interval < x < 1 which parametrizes the 
hierarchically organized matrix Qab, reflecting the ultramet- 
ric structure of the low-energy phase-space fTPi. In general, 
q{x) exhibits a plateau at large and small x, q{x > Xc) = Qc, 

q{x < X,n) = Qm- 

Jumps in the equilibrium configuration as a function of h are 
closely related to chaos in a field. Equilibrium configurations 
in different fields have minimum overlap as soon as the differ- 
ence in fields significantly exceeds 1/\/N [301, which sets a 
typical scale for large shocks. Interestingly, the same scale is 
also suggested by a dynamical consideration of local stability. 
The distribution of the local field hi — J2j^i JijO'j + h, i.e., the 
energy cost to flip only spin i, displays a linear pseudogap p2) , 
marginally satisfying the minimal requirement for metastabil- 
ity. The smallest local field thus scales as 1/^/N, setting the 
scale for Sh required to trigger a (dynamical) avalanche, as ar- 
gued and confirmed numerically in 126) . Accordingly, in both 
cases, the average size of an avalanche (total magnetization 
change) should scale as AAI ^ Nx6h ^ VN, where x is 
the average susceptibility. This is confirmed below. 
Although we now outline the main steps of the calculation on 
the SK model, the technique immediately extends to any mean- 
field system described by replica-symmetry breaking saddle 
points. Details will be presented in |33|. The probability for 
a shock in the interval [h, h + 6h/y/N] is proportional to 5h 
if (5/i <C 1; its fingerprint are non-analyticities in the moments 

of magnetization differences, [M{h) — M{h + Sh/y/N)]'' ^ 
N''/^\Sh\. Calculating the prefactor of |(5/i| for all k allows us 
to infer the avalanche-size density per unit field for Am > 0: 
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where we have introduced the suitably rescaled magnetization 
m ~ AI/N^/^, which jumps by Am ~ 0(1) in typical shocks. 
To calculate correlators of magnetization in different fields, 

Af (/ii) . . . M{hk) - (-l)'a^, . . . dh,F{hi) . . . F{hk), (3) 



we consider the generating function of a = 1. ..., n replica 
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Organizing the n replica into k groups subject to the same field 
^1=1,. ..,fc = h + hi/\/N, with ^^ ha = 0, and analyzing the 
cumulant expansion of the potential VF[{/ii}], the fc-point cor- 
relator ([3]) can be extracted in the limit n — > 0. Expanding 
A{Q, {ha}) to second order in hi, the potential is evaluated at 
the saddle point where Qab assumes Parisi's equilibrium so- 
lution qjiiy). However, due to the explicit breaking of replica 
symmetry by the external fields ha, a sum over inequivalent 
saddle points differing by replica permutations of Qab has to 
be performed. Generalizing techniques introduced in [34 1, we 
find a compact integral representation for the fc'th cumulant 
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where 0(x, y) solves the differential equation 
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In order to unambiguously identify shocks we need to take the 
Umit 7V"i/2 > T ^ 0. It is known |jT8]|20| that the non- 
analyticities oc | ft,,; | in the cumulants are obtained by an expan- 
sion of the diffusion-type equation dSl) to first order in the last, 
non-hnear term. For k > 2, the result encapsulates the full 
statistical information about jumps p3). 
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where /11.2 = /i + hi^2/\~N, and /112 — /ii— /12 > and a 

density (per unit of 5h) of jumps of size Am > 0, cf. Eq. ilV' 
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The weight t'^(g) = \iniT^o[Tdqjj^/ dx]^^ can be interpreted as 
the probability density, per unit energy, of finding a metastable 
state at overlap within [q, q + dq] with energy close to the 
ground state |31|. The density of shocks receives contribu- 
tions from the largest (q '^ qc{T — 0) — 1) to the small- 

2/3 

est overlaps qm{h) ~ h . Jumps in overlap of order 0(1) 

are indeed expected due to field chaos | [30) . A useful check of 
Eq. (J7]i is provided by the average magnetization jump which 
turns out to equal the thermodynamic (field cooled) susceptibil- 
ity, J pj^{Am)AmdAm — limT^oT^^ L dx{qc — (/(x)) ~ 
Xfc{T — 0). This is expected since the intra-state (zero-field 
cooled) susceptibility vanishes as T ^ 0, the susceptibility re- 
sponse being entirely due to interstate transitions. 
The formula Q has a very natural interpretation. If we take 
hi2 <C 1 in (rob we only need to consider the possibility that 



'it contains a piece S{q — qm)xm/T when ^(x) exhibits a plateau at x < 
I (if /i 7^ 0), hence the notation q^ in the integral. The integral measure can 



also be written as JT " ' d{x/T). 



the ground state and the lowest-lying metastable state cross as 
we tune h from hi to /i2, corrections being of order 0{hf2)- 
The disorder-averaged density of states of this two-level sys- 
tem is described by vj^{q)dq dE. The two states differ in N^ = 
N{l — q)/2 flipped spins. In the SK model the magnetization is 
uncorrected with the energy, and one thus expects the magne- 
tization difference between the states to be a Gaussian variable 
of zero mean and variance {Am?)q — ANn/N — 2(1 — q). If 
Am > 0, a jump at equilibrium occurs once hi2 = E/ Am. 
For the shock probability per unit h one thus expects 
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reproducing precisely Eq. (jTli. The above result (JTl) is generally 
valid for models described by RSB. It thus applies to p-spin 
models, where there is only one step of RSB, qj^^^ix) = qo + 
(91 ^ 'Zo)^(x — Xi). The avalanche distribution then simplifies 
with J dq v{q} — >■ Xi J dqd{q — qo), xi = Xi/T, into the form 



pi?>^'(Am) = XiAm^ 

^ V47r(g: 
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One verifies that its second moment agrees with Ref. p8| . The 
distribution (|9]) is non-critical, peaking around a typical size 
Am ^ 2y/qi — qQ, with p(Am) ^ Am at small Am (similar 
to one of the lower curves in Fig. [T]). The case of SK with full 
replica-symmetry breaking is much richer, as there is a T = 
limit function q{x). The weight with which events at overlap 
distance I — q contribute is a power-law |35 1, 



i^(g|l»l-g»T2)=C(l-9)-^/2 



(10) 



with C = 0.32047 ||36|. This holds independently of the exter- 
nal field h, and of additional random-field disorder [37]. From 
( 10 1 and (JTl) it leads to a robust scale-invariant jump density: 
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with T = 1. The universal exponent r = 1 for jump sizes 
7V^i/2 ^ Am, ^ 1 results from superimposed contributions 
from all overlaps, i.e. all scales, illustrated in Fig.fT] The cutoff 
function for larger jumps Am > 1 depends on the applied field. 
In zero field, ^(x) is linear at x <C 1. The resulting density 
i^(0) = 1.34523 at g = leads to the asymptotics 
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with t' = 1. Plots at intermediate Am = 0(1) are shown in 
Fig.fT] using approximations to q{x). A small field produces a 

plateau at (7min(^) ~ 1.0 x ft and while ( |ll| remains un- 
changed, the asymptotics ( 12 1 for Am S> Am/, ^ ft^^/"^ now 



decays with t' — —1, as for the one-step RSB case, replacing 
in^xi^Xh- iyiO)qmin(h) and qi -^ qniin(h)- 
Accepting Eq. dSll to represent the joint distribution of q and 
Am, we can integrate it over Am instead of q, which gives 
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Fig. 1: Power-law distribution of jumps for the SK model, from con- 
tributions from all overlap 1 — q. The curves in the lower part show 
the one-step like contributions for {1 — q) — 2~^ , fc = 1, ..., 12. The 
three nearly coinciding lines on the top show p(Am) from Eq. iItI, 
for h = 0,0.25 and 0.5, respectively. We use the approximation 
x(g) = {aq + hq^)l^\ - q with a = 1.28 and h = -0.64, and a 

2/3 ^ " " ' 

sharp lower cutoff at gniin(/i) ~ l-O/i |35,38,40|. The increase of 
h decreases the cutoff at large Am, while the avalanche distribution 
for Am ^ 1 is a universal power law. 



the probability distribution ViN^dN^ to flip A'a spins when 
increasing the magnetic field by 6h\ 



2? iVfl = 
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with p = 1. A very similar density of avalanches with the 
same exponents t = p = 1 was observed in the T = hys- 
teresis curve of 1 26] p In both cases, the number of spin flips 
scales as Nn ^ N"', with cr = 1, whereas the magnetization 
changes as Am ^ N^, with f3 = 1/2. This coincidence be- 
tween equilibrium and driven dynamics is presumably related 
to the marginality of the spin glass. It may also be due to the 
fact that the system is in a mean field limit, which, in the case 
of elastic manifolds at iV = 1, indeed gives the same exponents 
pT) . Similar coincidences were reported in other models f39l. 
At T = 0{N^^^^) > 0, non-analyticities in the moments are 
rounded |[T9]|28). An exact calculation for the second moment 
in the expansion of Eq. (|5]l shows |33| that the shock-related 
cusp ~ |/i| turns into an analytic crossover function f{x) of the 
scaling variable x = hi2/T, with J{x) -^ \x\ for large x. 
We calculate the distribution of shock sizes as an average over 
samples for a fixed small interval 5h ^ l,a priori not the same 
as varying ft. in a given sample. Such kind of self-averaging is 
expected for finite-dimensional glasses, whereas in the mean- 
field case it may be more problematic. It could be tested by 
computing correlations between subsequent shocks, extending 
our calculation to (5 /i = 0(1), or by analyzing the joint distri- 
bution of energies and magnetizations of low -lying states |40|. 
Finite-dimensional spin glasses - What aspects of our mean- 
field analysis are expected to survive in finite d-dimensional 

^A comparison of the numerical prefactor is unfortunately not possible, be- 
cause the avalanche density in I26j has been normalized by a factor ^ log(7V). 



spin glasses? We argue that independently of whether replica- 
symmetry breaking pTj or the droplet picture fTT) describes 
the glass state, at low fields the distribution of equilibrium 
avalanches is expected to be a power law. Indeed, let us assume 
that the dominant low-energy excitations are droplet-like spin 
clusters that flip simultaneously. These droplets are clusters 
that cannot be decomposed into a set of independent smaller 
excitations with lower energies. For droplets of typical linear 
size L we assume a typical energy cost L^ and a non-vanishing 
density of states (per unit volume) down to i^ = 0: vl {E = 
0)dE = vqL^'^^ L^^ dE with a constant vq independent of L. 
Empirically 6 is always very small, and df is the (possibly frac- 
tal) dimension of the droplets. We assume the total magnetiza- 
tion of droplets of size L to be uncorrected with the energy, and 
distributed a^PL(AAf) = L-'*»?/'A/(AA:^/i'^")- In a vanish- 
ing field, low energy droplets are believed to exist at all length 



scales, while recent numerical results 1 41 1 point towards the ab- 
sence of a thermodynamic glass phase in a finite field h. This 
implies a finite lengthscale Lh ^ 1/h beyond which droplets 
are suppressed. We make the standard assumption that droplets 
at scale L are uncorrelated with droplets at scales > 2L. With 
a reasoning analogous to the one leading to Eq. dSll, we expect 
a power-law density of avalanche sizes AM (per unit volume 
and unit field, with 6h — > 0): 
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(14) 
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with exponent r = ^^^ and a cut-off AM — .1^^ 
merical investigation of avalanches at small fields couH yield 
insight into the various exponents entering ( 14 1. Further- 
more, experimental measurements of power-law Barkhausen 
noise in spin glasses (e.g., by directly monitoring magnetiza- 
tion bursts 1 8 , 42 1) could provide complementary insight to ear- 
lier investigations of equilibrium noise |43 1. 
Elastic manifolds - The above calculation directly applies to the 
A^-component elastic manifold with coordinate u{x) of internal 
dimension d in a random potential, in presence of a harmonic 



well of curvature m , which forces (u) 



It has energy 
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As V is increased at T = along a straight line, i.e. Vi = vSn, 
the minimum-energy configuration jumps and static avalanches 
occur, of size S — J 6ui{x). We study models with Gaussian 
bare disorder V{u,x)V{0,0) ~ 6'^{x)Ro{u) with correlator 
Raiu) = NB{u^/N) and B'{z) = -(1 + ^)-'^, 7 > 0, in 
the large- A^ limit. Notations are as in [201, except here the 
Parisi variable is denoted x = Tx (u — Tu there). In |20| the 



^^We note that the numerical study (44( found that most likely neither of the 
assumptions dm = d{/2 = d/2, made by Fisher and Huse 1 17], holds. 

"^The mean-field case t = 1 is formally recovered replacing L'' — > (1 — q), 
dm/d — >• 1/2 and (df + 9)/d ^1/2. The latter reflects the typical gap at 
distance 1 - g. A, = (1 - q)^^^ jioi. 
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second moment of the renormalized disorder correlator, R{v) 
was computed. Here we are interested in the shock-size density 
p{S) — poP{S), the total shock density po and the normalized 
size distribution (/ dSP{S) — 1). We define its moments as 
(SP) = JdSSPP{S). The dictionary is as follows: h -^ vi. 
Am -^ m^S, g(x) -^ m^L'^G{x) with G(x) = G(/c = 0,x = 
Tx), where (u°^,u^) = Gab{k)- This gives the shock density 
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Two exact relations hold in all cases: 

f dSSpiS) = pa{S) = L'^fl - - 
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The first one is the total susceptibility 9^ J ui (x) = L'^ minus 
the intra-state susceptibility. The factor (1 — t-) thus gives 
the fraction of motion which occurs in jumps, which vanishes 



at m > m^ = mc(m°). Here trie 
Larkin mass, defined as m^ 



= mc(m) is the running 
[<T](?i+) in the notations 
of L20J . Eq. (1 8 1 extends the relation obtained in [11] between 
size moments and the cusp of the force correlator to the case of 
a finite fraction of motion in shocks. The size of the cusp is the 
same as in | ,20J . We define the large-size cutoff scale Sm via. 



G(x+) - 
For d < A and m < m" 
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= (Md)!/^, where e 
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A - d and 
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the T — Q^ saddle point equations admit a 



RSB solution |20 22j. We now discuss various cases, depend 



_ 2+7(d-2). 
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ing on the energy exponent 6 = 

(i) one-step RSB: it occurs for 9 < 0, i.e. d < 2 and 7 > ji 
The shock-size distribution depends on the single scale S^'- 
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with f{x) = a;''(x-2-e _ ^yfz ^ j^^ dy{y~^-y - l)i/2. For 
m <C rric, the size distribution becomes P{S) ~ ^P(^) with 
typical size Sm R 
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, avalanche-size exponent r 
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(23) 
*. One has (S) = ^/^^S'^, 



where r(a, z) — J 
and po = L'^/{S) from (17 1, consistent with (22 1 at small m. 
Note that for S ^ Sm, the first term, hence the one-step form, 
dominates in (|23|. There is however a distinct small-size cutoff 
scale, Sc = {-^)^^ Sm such that P{S) is a pure power law, 
P{S) - S-" for Sc < 5 < Sm- Since J^ dsp{s) ^ 1, the 
region S ^ Sm contains all the weight, consistent with r < 1 
(9 < 2). For sizes S ^ Sc , the probabiUty P{S) vanishes as 
P{S) ^ S with a peak around Sc- 



Interestingly, the droplet argument ( 14 1 can be adapted to the 
interface, i.e. A^ = 1. The correspondence AA/ — > S and 
standard interface scaling implies d^ ^^ d + C, and df —> d. 
Together with 6* == d - 2 + 2C it yields t = r^; = 2 - 2/(rf + C) 
and provides, for static avalanches, a basis for the conjecture 
made previously at depinning, i.e. out of equilibrium |45|. By 
contrast, the above large-A^ limit gives r = 2 — 2/(| + (). 
In d = 4, this gives t = 1, which is different from the usual 
mean field exponent t — 3/2 at A^ = 1 lllJ, and expected 
to hold at finite N. There are indications that this is due to a 
non-commutativity of the limits N.L -^ 00, also reflected by 
the unusual L-dependence of the maximal avalanche size 5m- 

Decaying Burgers - We now consider the decaying Burgers ve- 
locity field u(r. t) in dimension N, satisfying 



with Gaussian, power-law correlated, initial condition 

u,(r,t = 0)u,(r',i = 0) - -d,d,Roir - r') ^\r-r'\ 



-s'/i 



(20) 



(24) 

-27^ 
(25) 



Hence (5) = ^Sm which yields po from (17 1 



Here 



m^ ) depends on the details of the one- 
step solution. In the critical limit m ^ rric, for d > (d — 
is treated below), one finds 5m 

mc = [ — ''j^~ ] , and a roughness exponent C = (2 — d)/2 

(defined by u ^ x''). 

(ii) continuous RSB: it occurs for 9 > 0, with mc(m) = m^. 



ylG(x) = 



8 1 2 



Ro{r) ~ NB{r'^/N). From the Cole-Hopf transformation, the 
velocity at time t = 1/m^ is obtained from the d — version 
of the model ( 15 1 with v = r and T = 2i/ as tui{r, t) = Vi — 
{ui)-u- In the large-dimension limit iV — > 00 and for LR corre- 
lations < 7 < 1, i.e. < = 2^^ < 2, the above results for 
' 1+7 ' 

mj7^(mL)'^/^m^'^^'' with the manifold immediately apply, setting d = 0. In the inviscid 

limit, v ^f Q, the velocity field develops discontinuities along 
codimension-one manifolds, i.e. shocks, for t > tc = (m°)~^. 
Consider a line, say r^ — riSn, and the velocity field along 
this line i.e. ui(r) jumps by Au — Ui(7-j^) — Ui(rj^). From 
the identification Au = m^5, the shock-size density along this 



yl\i+2/e 

X 



and G(x) = G(xm) for x < Xm with Xc = Am* 

've \ g / 



Xm < X < Xc 

(21) 
Xm = Am", 



line is given by Eqs. ( 16 1 and (21 1 as p(Au) = m p(5) with 



1 20 



. The total shock density is 



m 



'■Ld/2 



PO 



2^ f-i^ 



(22) 



m = 1/t. At large time t ^ tc, the size probability takes the 
form P(Au) = (Auf)^^|5(Au/Aut) where p{s) is given by 
( [23] ) and the shock-size exponent is r = 1 — 7. The typical 
shock size is Au* = m^5m ^ i^^+2, with ( = 2/(1 + 7), 
consistent with the standard asymptotic scaling of the decay- 
ing velocity field: u{r,t) '" '"^> t^^+2u(f — rt~''^^) 



The 



p-5 
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total shock density po given by ( 22 1 vanishes for t < t^, ex- 
hibits a maximum near tc, then decays as ^ i^^/^, as the 
shock separation grows as t^/"^ from usual scaling. For shorter- 
ranged initial correlations, 7 > 1, the solution is one-step RSB 
and the reduced size distribution is p{s) 



:(1 



*c ' 



and (i) Au^ 



t 



-1/2, 



\se'''' /* with 
1 - i)i/2 for 



the 7 = 1 LR class; (ii) Aut ^ t~'^/^{^\nt)-^/^ ?A\aigtt, 
for the short-range class, very similar to the Kida result |24| for 
A^== 1. 

Conclusion - Systems whose thermodynamics is described 
by full RSB exhibit a power-law distribution of equilibrium- 
avalanche sizes, which can be traced back to their marginal sta- 
bility. Even though dynamic avalanches are different from our 
static analysis, the exponents turn out to be the same t = p = 1 
in the SK model, and in both cases the scale-free response is 
a consequence of criticality and marginal stability |26) . We 
expect a similar critical response upon slow changes of sys- 
tem parameters in many other systems with full RSB. This 
is of interest for optimization problems on dilute graphs such 
as minimal vertex cover |46|, coloring or Potts glass pT) , k- 
satisfiability |48| around the satisfiability threshold, and even 
in the whole UNSAT region at large k. Likewise, in models 
of complex economic systems one expects a power-law dis- 



tributed market response to changes in prices and stocks |49|. 
Finally, avalanches are expected in electron glasses with un- 
screened 1/r interactions. A stability argument shows that the 
number of rearrangements upon adding a new electron at T = 
diverges with system size at least as L'^^^, presumably with 
a wide distribution of dynamic responses. Since mean field 
yields a full RSB phase [37|, we speculate that static avalanches 
are power-law distributed as well. 

We thank S. Franz and M. B. Weissman for useful discussions. 
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REFERENCES 



[1] J. R Sethna, K. A. Dahmen, C. R. Myers, Nature 410, 242 

(2001). 
[2] T. Emig, R Claudin, J.R Bouchaud, Europhys. Lett. 50 594 

(2000) 
[3] D.S. Fisher, Rhys. Rep. 301, 1 13 (1998). 
[4] D. Bonamy, S. Santucci, L. Ponson, Rhys. Rev. Lett. 101, 

045501 (2008). L. Ronson, Rhys. Rev. Lett. 103, 055501 

(2009). 
[5] R Le Doussal, K.J. Wiese, S. Moulinet and E. RoUey, ERL 87 

(2009)56001. 
[6] E. Altshuler et al.,Rhys. Rev. B 70, 140505(R) (2004). 
[7] V. Repain et al., Europhys. Lett. 68, 460 (2004). 
[8] J. S. Urbach, R. C. Madison and J. T. Markert, Rhys. Rev. Lett. 

75, 276 (1995). D.-H. Kim, S.-B. Choe, and S.-C. Shin, Rhys. 

Rev. Lett. 90, 087203 (2003). 
[9] D. Monroe et al., Rhys. Rev. Lett. 59, 1148 (1987). M. Ben- 

Chorin, Z. Ovadyahu and M. RoUak, Rhys. Rev. B 48, 15025 

(1993). A. Vaknin, Z. Ovadyahu, and M. Pollak, Rhys. Rev. 

Lett. 81, 669 (1998). T. Grenet, J. Delahaye, M. Sabra and F 

Gay, Eur. Rhys. J. B 56, 183 (2007). 
[10] R Bak, C. Tang, and K. Wiesenfeld, Rhys. Rev. Lett. 59, 381 

(1987), D. Dhai-, Physica A 263, 4 (1999). 



[11] R Le Doussal, K.J. Wiese, Rhys. Rev. E 79, 051106 (2009). 
[12] R Le Doussal, A. A. Middleton and K.J. Wiese Rhys. Rev. E 79 

(2009) 050101 (R). 
[13] A. Rosso, R Le Doussal and K.J. Wiese, Rhys. Rev. B 80, 

144204 (2009). 
[14] J. Sethna et al., Rhys. Rev. Lett. 70, 3347 (1993). O. Rerkovic, 

K. Dahmen, and J. Sethna, Rhys. Rev. Lett. 75, 4528 (1995). 

K. Dahmen and J.R Sethna, Rhys. Rev. B 53 14872 (1996). 
[15] E. Vives et al. Rhys. Rev. E 52, R5, (1995). 
[16] C. De Dominicis and I. Kondor, Rhys. Rev. B, 27, 606 (1983). 
[17] D.S. Fisher and D. A. Huse, Rhys. Rev. B 38, 386, (1988). 
[18] L. Balents, J. R Bouchaud and M. Mezard, J. Rhys. I 6, 1007 

(1996). 
[19] R Le Doussal, K.J. Wiese, Rhys. Rev. B 68, 174202 (2003). 
[20] R Le Doussal, M. MuUer, and K. J. Wiese, Rhys. Rev. B 77, 

064203 (2008). 
[21] R Le Doussal, Annals of Rhysics 325 49 (2010). 
[22] M. Mezard and G. Rarisi, J. Rhys. I (France) 1, 809 (1991). 
[23] J. Bee and K. Khanin Rhys. Rep. 447, 1 (2007). 
[24] S. Kida, J. Fluid Mech. 93, 337 (1979). 
[25] J-R Bouchaud, and M. Mezard, J. Rhys. A 30, 7997 (1997). 
[26] F. Razmandi, G. Zarand, G. Zimanyi, Rhys. Rev. Lett. 83, 1034 

(1999). 
[27] A. R Young, A. J. Bray, and M. A. Moore, J. Rhys. C 17, L149 

(1984). 
[28] H. Yoshino and T. Rizzo, Rhys. Rev. B 77, 104429 (2008) 
[29] R Krzakala and O. C. Martin, Eur. Rhys. J. B 28, 199 (2002). 
[30] S. Franz and M. Ney-Nifle, J. Rhys. A 28, 2499 (1995). 
[31] M. Mezard, G. Rarisi and M.A. Virasoro, Spin-Glass Theory 

and Beyond, World Scientific, Singapore (1987). 
[32] D. J. Thouless, R W. Anderson, and R. G. Ralmer, Rhilos. Mag. 

35, 593 (1977). 
[33] R Le Doussal, M. Miiller, and K. J. Wiese, in preparation. 
[34] J. R Bouchaud, M. Mezard, and G. Rarisi, Rhys. Rev. E 52, 

3656 (1995). 
[35] G. Rarisi and G. Toulouse J. Physique (Paris), Lettres, 41, 

L361(1980). J. Vannimenus, G. Toulouse, G. Rarisi, J. de 

Physiquel, 42, 565(1981). 
[36] S. Pankov, Rhys. Rev. Lett. 96, 197204 (2006). 
[37] M. Miiller and S. Pankov, Rhys. Rev. B 75, 144201 (2007). M. 

Miiller and L. B. loffe, Rhys. Rev. Lett. 93, 256403 (2004). 
[38] R. Oppermann and M. J. Schmidt, Rhys. Rev. E 78, 061124 

(2008). 
[39] Y Liu and K. A. Dahmen Rhys. Rev. E 79, 061 124 (2009). 
[40] S. Franz and G. Rarisi, Eur. Rhys. J. B 18, 485 (2000). 
[41] A. R Young and H. G. Katzgraber, Rhys. Rev. Lett. 93, 207203 

(2004). T. Jorg, H. G. Katzgraber, and F Krzakala Rhys. Rev. 

Lett. 100, 197202 (2008). 
[42] K. Komatsu et al., AIR Conf. Proc, 1129, 153-156 (2009). 
[43] M. B. Weissman, Rev. Mod. Rhys. 65, 829 (1993). 
[44] J. Lamarcq, J. -P. Bouchaud, and O. C. Martin, Rhys. Rev. B 

68, 012404 (2003). 
[45] O. Narayan and D.S. Fisher, Rhys. Rev. B 48 (1993) 7030. 
[46] A. K. Hartmann and M. Weigt, J. Rhys. A 36, 1 1069 (2003), R 

Zhang, Y. Zeng, and H. Zhou, Rhys. Rev. E 80, 021 122 (2009). 
[47] R Krzakala and L. Zdeborova, ERL 81, 57005 (2008). 
[48] A. Montanari, G. Rarisi, and F. Ricci-Tersenghi, J. Rhys. A 37, 

2073 (2004). 
[49] J. R Bouchaud, Rhys. World, April 2009, pp 28-32. 



p-6 



